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ABSTRACT 

Heavy Quark Effective Theory (HQET) is a new approach to 
QCD problems involving a heavy quark. In the leading approx- 
imation, the heavy quark is considered as a static source of the 
gluon field; 1/ra corrections can be systematically included in the 
perturbation theory. New symmetry properties not apparent in 
QCD appear in HQET. They are used, in particular, to obtain re- 
lations among heavy hadron form factors. HQET also simplifies 
lattice simulation and sum rules analysis of heavy hadrons. 

Part 1 contains discussion of the effective lagrangian, mesons, 
baryons, and renormalization. Part 2 will contain 1/m correc- 
tions, nonleptonic decays, and interaction with soft pions. 
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1. Effective Lagrangian 



Recently an interesting new approach to QCD problems involving a heavy 
quark was proposed, namely the Heavy Quark Effective Theory (HQET). In 
the leading approximation, the heavy quark is considered as a static source 
of the gluon field; 1/m corrections can be systematically included in the 
perturbation theory. This simplification is similar to considering a hydrogen 
atom instead of a positronium. New symmetry properties not apparent in 
QCD appear in HQET. They are used, in particular, to obtain relations 
among heavy hadron form factors. However, in QCD even such a simplified 
problem is unsolvable. Approximate methods such as lattice simulation or 
sum rules are necessary to obtain quantitative results. Here again HQET 
allows to proceed much further than QCD. 

There are several good reviews of HQET 0HII to which we address the 
reader for an additional information. Here we widely use the properties of 
currents' correlators to obtain general results. This approach is inspired by 
sum rules, though we shall not consider details of sum rules calculations. We 
shall start from a very simple though approximate picture in the Sections |l| ^: 
some complications are discussed later. 

Let's start from the QCD Lagrangian 

L = Q(iD - m)Q + qiDq - \g%G% + ■■■ (1.1) 

where Q is the heavy quark field, q are light quark fields (their masses are not 
written down for simplicity), G" is the gluon field strength, and dots mean 



3 



gauge fixing and ghost terms. It is well known that the free heavy quark 
Lagrangian Q(id — m)Q gives the dependence of the energy on the momen- 
tum e = y 'm 2 + p 2 . We shall consider problems with a single heavy quark 
approximately at rest, and all characteristic momenta \p\ <C m. Then we 
can simplify the dispersion law to e = m. It corresponds to the Lagrangian 
Q(z7o(9o — m)Q. In such problems it is convenient to measure all energies 
relative to the level m. This means that instead of the true energy e we shall 
use the effective energy e = e — m. Then the heavy quark energy e — inde- 
pendently on the momentum. The free Lagrangian giving such a dispersion 
law is QijodoQ. The spin of the heavy quark at rest can be described by a 
2-component spinor Q (we can also consider it as a 4-component spinor with 
the vanishing lower components: 70Q = Q)- Reintroducing the interaction 
with the gluon field by requirement of the gauge invariance, we arrive at the 
HQET Lagrangian [| 

L = Q + iD Q + qiDq-jG° v G% + --- (1.2) 

The static quark held Q contains only annihilation operators. There are no 
heavy antiquarks in the theory, because processes of heavy quark-antiquark 
pair production are suppressed by 1/m. The heavy antiquark (if present) is 



described by a separate held. The held theory (L2) is not Lorentz-invariant 



because the heavy quark defines a selected frame — its rest frame. 



The Lagrangian (1.2) gives the static quark propagator 



S(P) = ^-^> S(x)=S(x )6(x), S(t) = -«(t). (1.3) 
Po + «0 

In the momentum space it depends only on p but not on p because we have 
neglected the kinetic energy. Therefore in the coordinate space the static 



quark does not move. The unit 2x2 matrix is assumed in the propagator (1.3) 



It is often convenient to use it as a 4 x 4 matrix; in such a case the projector 
1+ 2° excluding the lower components is assumed. The static quark interacts 
only with A ; the vertex is igS 0lJ ,t a . 

One can watch how expressions for QCD diagrams tend to the corre- 
sponding HQET expressions in the limit m — > 00 [||. The QCD heavy quark 
propagator is 

of, P + m m(l + 70) +p 1 + 7Q , n ( P \ n 

S (P) = ^ 2= 0~i~2 = 0~ + ~ • L4 

pz _ m z zmpo + p z 2po \ m ) 

A vertex igj^t a sandwiched between two projectors 1+ 2° may be replaced 
by ig5f)^t a (one may insert the projectors at external heavy quark legs too). 
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Therefore any tree QCD diagram equals the corresponding HQET one up 
to 0(p/m) terms. In loops, momenta can be arbitrarily large, and the re- 



lation (1.4) can break. But regions of large loop momenta are excluded by 
the renormalization in both theories, and for convergent integrals one may 
use (see Sec. §. 



The Lagrangian ([L^) can be rewritten in covariant notations: 

L = Qiv^D^Q + --- (1.5) 

where the static quark field Q is a 4-component spinor obeying the relation 
vQ = Q and w M is the quark velocity. The true total momentum p^ is related 
to the effective one p^ by 

Pn = TnVft + p^, \pn\ < to. (1.6) 

The static quark propagator is 

and the vertex is igv^t a . In the limit m — > oo the heavy quark can't change 
its velocity in any processes with bounded momenta p^. Therefore there 
exists the velocity superselection rule Q: heavy quarks with each velocity 
Vp_ can be treated separately and described by a separate field Q v . If we are 
interested in a transition of a heavy hadron with the velocity v\ into a heavy 
hadron with the velocity vi , we can use the Lagrangian 

L = 5^Q i w iM I> M Qi + "- (1.8) 

i 

where Qi is the static quark field with the velocity x>i (the quark Q\ is present 
in the initial hadron and Q 2 — in the final one). These quarks have different 



propagators (1.7) and vertices. They can be of the same or different flavour; 
it doesn't matter because they can't transform into each other except by an 
external current with an unbounded momentum transfer (of order to). It is 
even possible to write a Lorentz-invariant Lagrangian 

L = J ^Q v iv^D^Q v + --- (1.9) 

describing static quarks with all possible velocities at ones. But in any specific 
problem only several heavy quarks with several velocities are involved; all 
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fields Q v except few ones are in the vacuum state and are irrelevant, and 
finite sums ( |l.§| ) are sufficient. 

There is an ambiguity what quark mass m should be used in ( jl.6| ) [Q . In 

general the HQET Lagrangian is Qiiv^D^ — Sm)Q; the residual mass Sm is 
shifted when we change m. Physical quantities, of course, don't depend on 
this choice. The most convenient definition of the heavy quark mass m is one 
that gives 5m = 0. It corresponds to the pole of the quark propagator at 
VfiPn = 0, or p 2 — m 2 in QCD. This pole mass is gauge- invariant. There is 



also an ambiguity in the exact choice of w M in (1.6) fel. This reparametrization 



invariance relates coefficients of terms of different orders in 1/m expansion. 



Quantization of the theory (1.8) was discussed in 



The HQET Lagrangian (1.2) possesses the SU(2) spin symmetry [lCj ]. 
The heavy quark spin does not interact with gluon field in the limit m — > oo 
because its chromomagnetic moment vanishes. If there are rih heavy quark 
flavours with the same velocity, there is the SU(2rih) spin-flavour symme- 
try. For example, in the problem of transition from a heavy hadron with 
the velocity v\ to a different flavour he avy hadron with the velocity i>2, at 



equal velocities v\ — v 2 the Lagrangian (|l.8|) has the SU(4) spin-flavour sym- 
metry which relates all form factors to the form factor of a single hadron 
at zero momentum transfer (equal 1). At non-equal velocities, it has only 
the SU (2) x SU (2) spin symmetry relating form factors to each other. The 



Lagrangian ( |l.9|) has the symmetry SU(2rih)° c x 50(3, 1). 

Not only the orientation but also the magnitude of the heavy quark spin 
is irrelevant in HQET. This leads to a supersymmetry group called the su- 
perflavour symmetry |]llf . It allows one to predict properties of hadrons with 
a scalar or vector heavy quark appearing in supersymmetric extensions of the 
Standard Model, in technicolor models, and in some composite models. The 
scalar and vector static quark Lagrangians 

L = ip + iD !f + ■ ■ ■ , L = ViD V + --- (1.10) 

have the SU(rih) and SU^n^) spin-flavour symmetry. This idea can also be 
applied to baryons with two heavy quarks [ p^[ because they form a small size 
(of order l/ma s ) spin or 1 bound state antitriplet in color. 

HQET has great advantages over QCD in lattice simulation of heavy quark 
problems. Indeed, the applicability conditions of the lattice approximation 
to problems with light hadrons are that the lattice spacing is much less than 
the characteristic hadron size, and the total lattice length is much larger than 
this size. For simulation of QCD with a heavy quark, the lattice spacing must 
be much less than the heavy quark Compton wavelength 1/m. For b quark it 
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is impossible at present. The HQET Lagrangian does not involve the heavy 
quark mass to, and the applicability conditions of the lattice approximation 
are the same as for light hadrons |0. Relation of the lattice HQET to the 
continuum one was investigated in [ jl4]Jll||l(| . Simulation results can be found 
in 0. 

2. Mesons 

Due to the heavy quark spin symmetry, hadrons may be classified according 
to the light fields' angular momentum and parity j which are conserved 
quantum numbers. In other words, we can switch off the heavy quark spin 
using the superflavour symmetry, and then the hadron's momentum and par- 
ity will be j p . The Qq mesons are the QCD analog of the hydrogen atom. 
The ground-state (S-wave) meson has j p = the excited P-wave mesons 
have j p = 5 and | . When we switch the heavy quark spin on, each of 
these mesons becomes a degenerate doublet. Its components are transformed 
into each other by the heavy quark spin symmetry operations. We have the 
ground-state doublet CP, 1~, and the excited P-wave doublets + , 1 + , and 
1 + , 2 + . Splittings in these doublets (hyperfine splittings) are due to the heavy 
quark chromomagnetic moment interaction violating the spin symmetry, and 
are proportional to 1/to (Sect. 5). 

Form factors of the ground-state mesons in HQET were considered 
in |]l0| , p^|Jl9tl ; applications to semileptonic B decays were discussed in the 
review pc| and papers and to e + e~ annihilation — in J^2j2^]. Transition 
form factors to the P-wave mesons were considered in |24f25|. A general 
method of counting independent form factors applicable both to mesons and 
baryons was proposed in p^| , and an elegant explicit construction — in 27]. 
It was applied to ground state to arbitrary excited meson transitions in 2S ] . 
Two-point HQET sum rules were investigated in pg| , |30| , and three-point 
ones — in [^lj. Mesons in two-dimensional QCD with the large number of 
colors were considered in |}2| . Here we shall use correlators of currents with 
the quantum numbers of mesons in order to investigate properties of mesons 
in HQET. 

When the heavy quark is scalar, there is one bilinear heavy-light current 
without derivatives J s = Qfq (Qf is the heavy antiquark field). It has no 
definite parity; the currents J± = 2 7 ' J s have the parity P = ±1 because the 
P-conjugation acts as q — * joq. The current j+ has the quantum numbers of 
the ground-state j meson, and j_ — of the P-wave^ meson. Currents with 
the quantum numbers of mesons with higher j necessarily involve derivatives. 



7 



In the case of real-world spin i heavy quark, there are 4 bilinear currents 

without derivatives J — QTq. Indeed, because of 70Q = Q, the current with 
T = 70 reduces to T = 1; 7075 — to 75; a i — to ^\ — to s i3 klkj5- We are 
left with r = 75, 7 and 1, 775. The first pair with T anticommuting with 
70 has the quantum numbers of the ground-state - , 1~ doublet; the second 
pair with T commuting with 70 — of the P-wave + , 1 + doublet. 



x 




Figure 1: Correlator of two HQET heavy-light currents 

A correlator of any two currents containing the static quark field has the 
form (Fig. |l|) 

i<TJ 2 (x)j£(0)> = 5(x)U(xo), (2.1) 
1%) = J U^e^dt, IL(t) = J n(w)e- Mt ^. 

It obeys the dispersion representation 

00 00 
= f P( £ ) d£ + . n(t) = r p^e-^dcj + ■■■ (2.2) 

J £ - LU - lO J 


A subtraction polynomial in n(u) (denoted by dots) gives S(t) and its deriva- 
tives in II(i). We can analytically continue a correlator from the half-axis 
t > to imaginary t = —it. Then II(r) and p(u>) are related by the Laplace 
transform 

00 a+200 
n(r) = i J pi^e-^dcu, p{uu) = J n(r)e a,r dr, (2.3) 

a— ioo 

where a is to the right from all singularities of II(t). 
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The contribution of an intermediate state \h> with the energy e to II(£), 
n(w), p(w) is 

U h (t) = t<0\j 2 \h>iS(t)e- t7t <h\J+\0>, 

U h (u) = i<0\J 2 \h> I— -<fc|3+|0>, 

uj — e + iQ 

PH(u) = <0\n\h><h\Ji\Q>5{LO-e). (2.4) 

We remind the reader that the HQET energy e means the true energy minus 
the heavy quark mass. 

The correlator of two meson currents with the scalar heavy quark has the 
7-matrix structure (Fig. 

i<TJ s {x)JM> = 8{x)U s (x ), U s = A + B 7o . (2.5) 

For the currents with the definite parity P we have 

t<TJ P (x)J P (0)> =6{x)Tlp(x Q ) 1+ p ,Ilp = A + PB= ^ Tr(l + P 7o )II s . 

(2.6) 



Due to the linear relations (2. 1- 2. 3D , the same 7-matrix structures and re- 
lations between Tl s and lip hold in both the coordinate space and the mo- 
mentum one, and also for spectral densities. When calculating the correlator 
using the Operator Product Expansion (OPE), even-dimensional terms con- 
tain an odd number of 7-matrices along the light quark line and after all 
integrations contribute to B; odd-dimensional terms contain even number 
of 7-matrices and contribute to A. If we denote <0|J + |M, \ + > = f M i + u 
where u is the meson M wave function, then the meson's contribution to 
ps{w) summed over polarizations is /?. i + 8(u> — e) Y] uu, or the contribution 

to p+(u>) is f^, i_ + 6(u> — e). Similar formulae hold for i mesons. 

Now let's switch the heavy quark spin on. The correlator is (Fig. |l|) 

i<Tj 2 (x)jl(o)> ^ 8(x)u(x ), n = Trr 2 i^r 1 n s . (2.7) 

In II S , 70 may be replaced by P = ±1 for (anti-) commuting with 70, 
and IT becomes the scalar function lip: 

1 — 'Yn 

n = n P Trr 2 — — ri. (2.8) 
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The correlators of the currents Q"/§q and Qjiq with the quantum numbers of 
the ground-state 0~, 1~ mesons are equal to 2II + and 25ijH+. If we denote 

<0|Q7 5 g|Af,CT> = f M ,o-, <0|Q7g|M,l"> = f M ,i-e (where e is the 1~ 
meson's polarization vector) then the meson's contributions to the spectral 
densities are /| f0 _5(w — e -) ancl ^vImi-H^ — ^l-)- Therefore the spin 
symmetry requires that the mesons in 0~ and l - channels are degenerate 
(e - = el- = e)) and ./m,o- = Jm,i- = V^/^ i+- Similar formulae hold for 
P-wave + , 1 + mesons. 

In QCD, the meson constants are usually defined as <0| Q7 M 7sQ |M, 0~> = 
/mo-P(U> <0|Q7^<7|M, 1~> = mfj^i-e^, where the meson states are nor- 
malized in the relativistic way: <M,p'\M,p> = 2po5(p' — p). This nor- 
malization is senseless in HQET; we must use the non-relativistic nor- 
malization <M,p'\M,p> = S(p' — p) instead. Then the definitions read 
v / 2^<0|Q7o759l^,0-> = mf Mi0 -, V2m<0\Q^q\M, 1"> = mf Mil -e. Fi- 
nally we obtain the scaling law 



2/ 



M, 



1 + 



/A/,0- 



= f Mtl - = (2.9) 




Figure 2: Correlator of two HQET heavy-light currents and a heavy-heavy 
current 

To investigate hadron form factors in HQET, we^consider correlators of 
two currents J\ t2 containing the static quark fields Qi,2 with the velocities 
i»i,2 and the heavy-heavy velocity-changing current J (Fig. |J): 

oo oo 

t 2 <TJ 2 {x 2 )J{0)j+{x 1 )> = J dt 2 5(x 2 - v 2 t 2 ) J dtiSixi+vit^Kfah), 
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K(u 2 ,u 1 ) = J K{t 2 M)e W2t2+Wltl dt 2 dt u (2.10) 
K(t 2 , tl ) = j K{u», ^e-^-^ ^ ^. 

They obey the double dispersion representation 

s f p(e 2 ,si)de2dei 
K(u>2,ui) = / — 7 ttt 777 + ■ • ■ (2.11) 

K{t 2 ,t 1 ) = S{t 2 )S(t 1 ) j p{w 2 ,Lu 1 )e- lu)2t2 ~ luJltl dLu 2 duj l + --- 

Subtraction terms in K(u 2 ,u>i) (denoted by dots) are polynomial in u>i with 
coefficients that are arbitrary functions of ui 2 (given by single dispersion in- 
tegrals) plus vice versa. These terms give 8(ti) and its derivatives times 
arbitrary functions of t 2 plus vice versa in K (t 2 , t\). We can analytically con- 
tinue K(t 2 ,t\) from ti, 2 > to t\, 2 = —iri i2 . Then K(t 2 ,t\) and p(uj 2 ,uji) 
are related by the double Laplace transform 

K(T 2 ,n) = - J p{w 2 ,w 1 )e- ul2T2 - UJ ^dLJ 2 diJ 1 , (2.12) 

a+zoo a+ioo 

p(a, 2 , Wl ) = _I_ J dr 2 J d Tl Kfa,T 1 )e u > n+uin . 

a— too a— ioo 

The contribution of intermediate states |/ii,2> with the energies £1,2 to 
K(t 2 ,ti), K(u> 2 ,u>i), and p(u> 2 ,u>i) is 

K h2hl (t 2 ,h) - J 2 <0|J 2 |/ l2 >^(t 2 )e^ 72t2 </i 2 |J|/i 1 >z5(t 1 )e^* 1 <^ 1 |}+|0>, 

^ h2hl (w2,wi) - i 2 <0|J 2 |/i 2 > 1 , n <fe 2 |J|fei> ~ , . n < fe i|jil0>> 

^2 — £2 + Wl — £1 + IV 

Ph^uj^w-C) = <Q\J 2 \h 2 ><h 2 \J\h 1 ><h 1 \Jt\v>5{u 2 -e 2 )8{w 1 -ei), (2.13) 

where the sum over /ii ;2 polarizations is assumed. Let's introduce the spin 
wave functions ^1,2 of fti )2 . Then <0|Ji i2 |/ii,2> = f 1,2^1,2, <h 2 \J\h\> = 
ip 2 f 2 iipi, where /21 is a form factor matrix in the spin space. The spec- 
tral density of the correlator of the currents j^~ 2 V'i,2 with some specific 
polarizations Vi,2 is il> 2 p(uj 2 , u)i)ipi- The contribution of |/ii,2> to it is 

f2fl^ 2 f21^lS(uJ 2 - £ 2 )<5(wi - £1). 
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The correlator of two meson currents with the scalar heavy quark j s2 , J sl 
and the scalar heavy- heavy current J s — Q^Q S 2 has the 7-matrix structure 
(Fig-!) 

K s = A + Bivi + B 2 v 2 + Cv 2 vi. (2.14) 
For the currents Jp 2 , J Pl with the definite parities we have 

K P2Pl =A + P 1 B 1 + P2B2 + P1P2C = 4 \ ■ ^ l*l_L. 2.15 

1 + P1P2V1 ■ V 2 



Due to the linear relations (2.10-2.12), the same 7-matrix structures and 



relations between K s and Kp 2 p 1 hold in both the coordinate space and the 
momentum one, and also for spectral densities. When calculating the cor- 
relator using the OPE, even-dimensional terms contribute to B\, B2, and 
odd-dimensional ones — to A, C. 

It is convenient to use the "brick wall" frame in which v\ — — tf 2 is directed 
along z for counting form factors. Angular momentum projection onto z is 
conserved: j 2z = j\ z . The reflection in any plane containing z transforms 
the state \j,j z > to Pi 2 ^\j, — j z >. Therefore the amplitude for —ji x , —j 2z is 
equal to that for j\ z , j 2z up to a phase factor; the 0^0 transition is allowed 
only if the "naturalness" P{— 1) J is conserved ^(|. For example, | + — > i + , 
h — > I , and \ — > | transitions are described by one form factor each: 

1+ ~ 1 + 
<M,- \J S \M,- > = e(ch^)n 2 tti, 

1- ~ 1 + 

<M,~ \Js\M,- > = r 1/2 (ch^)u 275Ul , (2.16) 
3- ~ 1 + 

<M >2 l Js l M '2 > = T 3/2( chl ^) W lM"2 / i'Ul, 

where ch <p = V\ ■ V2 is the cosine of the Minkovskian angle between the world 
lines of the incoming and the outgoing heavy quark, and is the Rarita- 
Schwinger wave function of the spin | meson. Here we have slightly changed 
the notations as compared to p4| : in the original notations right-hand sides 



of ( 2.16 ) should contain 2ti/2 and V3t 3 / 2 - The contribution of | mesons 



to the spectral density p s (u)2,vi) is /2/i£(ch (^±^±4^(5(^2 - £2)^1 - £1), 
i. e. the contribution to p++(w2,^i) is /2/i£(ch ip)S(uj2 — £2)^1 — £i)- The 



spectral density of the correlator of the currents J s i, 2^1,2 with some specific 
polarizations ui 2 is U2ps(^>2> u>i)ui. The contribution of \ + mesons to it is 
/2Au 2 £(ch</3)ui<5(a->2 — £2)^(^1 — £1). Similar formulae hold for h mesons. 
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Now let's switch the heavy quark spin on. The correlator of j 2 = Q 2 T 2 q, 
h = qTiQi, and J = Q X TQ 2 is (Fig. |) 

K = TtT 2 ^-^T^-^T 1 K s . (2.17) 

In K s , W1.2 may be replaced by P\ 2 — ±1 for Ti.2 (anti-) commuting with 
W1.2, and K s becomes the scalar function Kp 2 p 1 : 

K = K P2 p 1 TrT 2 —^-T—£±T 1 . (2.18) 

Let's introduce the currents qM.Q, Ai = Tip where tp is the spin wave func- 
tion, i. e. M.q- = 7 5 = —75, Aii- = 7 p e M — e. The spectral density of 

their correlator is pp 2 p 1 Tr Ai 2 1 ~ 2 2 TAii; the mesons' contribution to it is 
7m 2 Jm 1 <M 2 1 J\Mi> 5(u 2 - e 2 )S(u}i - ei). Hence we obtain 

<M 2 |J|M 1 > = l^ c hip)TTM 2 ^-^T^^M 1 . (2.19) 

This formula expresses all form factors of transitions of a ground-state 0~, 
1~ meson to a ground-state 0~, 1~ meson under the action of any heavy- 
heavy current J = Q 2 TQi via one universal Isgur-Wise form factor £(chip). 
A similar formula with Ti/ 2 (chtp) holds for i + — > \ transitions; M. = 1, (375 
for + , 1 + mesons. 

A three-point correlator at <p = is expressed via the corresponding two- 
point one: 

K(t 2 ,h) = U(tx+t 2 ), 

K{ W ) = n(cJl) - nM , (2.20) 
Wi — u> 2 

The first two forms follow from each other by the Fourier transform. The 



third form is necessary and sufficient for the first one because of (2.2); it can 
be obtained by the double backward Laplace transform ( 2.12| ) of the single 



forward Laplace transform (2.3) of the first form at imaginary times. The 
third form can be also obtained from the second one by taking the double 
discontinuity: the discontinuity of the first term in u 2 is Tri5(uj 2 — u>i)IL(u>i), 
and the discontinuity of this expression in w\ is j(2tti) 2 p{ui)5(uj 2 — ^1); the 
second term contributes equally. 
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To prove the first form, let's consider any diagram for the two-point cor- 
relator in the coordinate space (for simplicity, with the scalar heavy quark). 
Vertices along the heavy quark line have the times to < ti < • • • < t n -i < tni 
and the integration in t^, ■ ■ ■ ,t n -i is performed. The integrand is an inte- 
gral over coordinates of all vertices not belonging to the heavy quark line. 
Now consider all diagrams for the three-point correlator obtained by insert- 
ing the heavy-heavy vertex with time t (and tp = 0) to all the possible places 
along the heavy quark line. These diagrams have the same integrand and 
the integration regions to < t\ < ■ ■ ■ < t m _i < t < t m < ■ ■ ■ < t n _i < t n 
(m = 1, . . . , n). These regions span the whole integration region of the orig- 
inal diagram. Therefore the sum of this set of three-point diagrams is equal 
to the two-point diagram. 

The second form can be easily proved in the momentum space in the exact 
analogy with the QED Ward identity using the relation iS(ui + uj')US(uJ2 + 

= l " ( Z\^ 2) (Fig- !)• In Particular, it implies K(cu,oj) = ^M. 
Comparing the mesons' contributions to p(w2,o-'i) and p(to), we see 

£(1) = 1 (2-21) 

for any i + meson. For non-diagonal ^ + — > i + transitions £(1) = because 
/?(il>2, i^i) = off the diagonal. The physical meaning of this is simple: when 
the current J replaces the old heavy quark by the new one with the same 
velocity and color, light fields don't notic e it. The formulae for the form 
factors at ip = equivalent to ( [2.19| , 2.21 ) were first proposed in the quark 



model framework |p3| . 

The variable ch ip is related to the momentum transfer q 2 (for simplicity 
in the case of equal heavy quark masses) by the formula q 2 — 2m 2 (1 — ch^j). 
Form factors are analytic functions of q 2 with the cut in the annihilation chan- 
nel from 4m 2 to +oo. Therefore the Isgur-Wise function £(chiy9) is an analytic 
function in the chip plane with the cut from —1 to — oo. Geometrically speak- 
ing, chip > 1 corresponds to Minkovskian angles between the incoming and 
outgoing heavy quark world lines (scattering or decay); ch^ = 1 means the 
straight world line — no transition at all; nothing special happens at chip < 1, 
only the angle becomes Euclidean; chi^ = — 1 is really a singular point where 
the quark returns along the same world line; chip < —1 corresponds again 
to Minkovskian angles only one of the world lines is directed to the past 
(annihilation). 

In the rest of this Section, we shall for simplicity live in the world with the 
scalar b quark decaying into the scalar c quark plus the scalar W boson. Their 
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Figure 3: Proof of the Ward identity. A digit 1 (2) near a heavy quark line 
means that its energy includes u\ {0J2) 
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masses can be adjusted in such a way as to give any desired ch ip. The quark 
decay matrix element is simply M — g where g is the coupling constant. 

Until now we discussed exclusive decays of the i + B meson. Inclu- 
sive decays can be also treated in HQET [ ^J3^| . The matrix element of 
the decay B — > X + W (where X is any hadronic state containing the 
c quark) has the structure M — gip 2 {X)u\. Its spin-averaged square is 

\M\ 2 = ^Tr±±^-0 2 (X)?A 2 (X). Let's sum over hadronic states X with the 
energy e: 

W(e,ch<p) = J2MX)^ 2 (X)S(£x -e)=a + bd 2 , 
x 

dY = dr w(e,chip)de, (2.22) 

w(e, ch if) = i Tr - W(e, ch(p) = a + b ch ip. 

The meson decay rate dT differs from the quark decay rate dTg by the struc- 
ture function w(e, chip) where £ is the energy of the hadronic state X (minus 
the c quark mass) in the v 2 rest frame. The quark-hadron duality tells us 
that the total meson decay rate is equal to the quark one. This is the Bjorken 
sum rule @|| 

OO 

J w{e,chip)de = 1. (2.23) 
o 

The second sum rule follows from the momentum conservation |}6]]. The 
initial ground-state meson has the energy e g in the v\ rest frame. In the v 2 
rest frame this corresponds to the energy e g ch ip plus an irrelevant momentum 
orthogonal to v 2 . Therefore the average energy of the hadronic system X must 
be 

OO 

J w(e, chip)ede = e g chp. (2.24) 
o 

Inclusive semileptonic B decays in HQET were also discussed in ]37[ . 

Now we shall explicitly write down some contribution to this sum rule. 
The spin-averaged matrix elements squared for the decays h + ~ ~ ' ' 



i , and i + -» | fl2T6|) 



2 2 ' 2 

are 



= ^Tvi±^i±^ = ^^±l, 

= ^Tr^ 75 ^ 5 = 9 ^ /2 ^1, (2.25) 
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Here we have used the Rarita-Schwinger density matrix. The decay „ — ► 5 
is S'-wave, hence |M| 2 is constant at ip — > 0. The decays ^ + — > 5 , ^ + — > | 
are P-wave, hence |M| 2 ~ tp 2 . The decays to the D-wave mesons | + , | + arc 
D-wave with |M| 2 ~ (p 4 , etc. Matr ix elements squared in the annihilation 
channel W — > differ from ( 2.25 ) by the absence of the factor i coming 



from the averaging over the initial meson spin states. The first decay W 



1 + 



is P-wave, the second one W 



is S'-wave, and the third one 



W - 
chip 



1-3- 
2 2 
-> -1. 



is £>-wave. This determines the threshold behavior of (2.25) at 



Therefore the Bjorken sum rule reads |24l 



^ 2 ch <g + 1 x - 2 chip - 1 x - 



2 (chv? + l) 2 (chp - 1) 

r 3/2 o 



• = 1, 



(2.26) 

where the sums are over resonances in the j p channels, and dots mean con- 
tributions of other channels (see |28[ ). The Burkardt sum rule ( [j.24[ ) has the 
similar form. The simplest consequence of (2.26) is that the decay rate to the 
ground-state \ meson is less than the total one: 



i{chip) < 



ch ip + 1 



(2.27) 



Of course, at chip^S> 1 most decays are inelastic, and £(ch ip) <C 1/y/chip. 

Let's consider the Bjorken sum rule ( [2.26 ) at small cht^ — 1, and expand 
it to the linear terms. The channels denoted by dots don't contribute be- 
cause they are at least D-wave. Higher resonances in the ^ + channel don't 
contribute because they have £(chy>) = 0(chip — 1). We are left with p3 



(2.28) 



This gives us the Bjorken bound £'(1 ) < — 4 (evident also from (2.27)). Sim- 
ilarly, the Burkardt sum rule (2.24) leads to the optical (Thomas- Reiche- 
Kuhn) sum rule pi] 



1 



4 ^ 



(ei/2 



)r 1 2 /2 (l) + 



-Y 
3 ^ 



( e 3/2 



£ g) T 3/2 



(1) = \e 9 - 



(2.29) 
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It can be used for obtaining bounds on £'(1) | pq |. 

It is also possible to establish the bound on the Isgur-Wise form factor 
at the cut the decay rate W — > BD is less than the total decay rate 



W — ► be. The meson decay rate for each flavour is given by (2.25) without 



the spin-averaging factor i; the quark decay rate is \M\ 2 = g N c where N c 
is the number of colors. If there are m light flavours for which £(chy>) is 
approximately the same, then 

ni\^{chip)\ 2 \chip+l\ < N c . (2.30) 

In general the left-hand side is the sum over light flavours. The factor ni 
was erroneously omitted in | ]39|; the phrase justifying this seems to have no 
sense. At |chy| ^> 1, the BD channel constitutes a small fraction of the 
total W — > be width, and |£(chy>)| <C 1/yj chip\. One could include also 
higher states' contributions in the left-hand side; this should be done with 
caution because of the possibility of double counting. The inequality ( 2.3C| ) 



is applicable only sufficiently far from the threshold, at | ch^ + 1| 3> 7r 2 «s- 
Near the threshold the Coulomb interaction between the heavy quark and 
antiquark is essential. The total decay width on the right-hand side is not 
equal to its free-quark value N c ; it contains high narrow resonances at the 
quarkonium levels. Moreover, the very concept of the Isgur-Wise form factor 
is inapplicable in this region. The HQET picture is based on the fact that 
heavy quarks move along straight world lines, but at velocities ~ Tta s they 
really rotate around each other. 



If the inequality (2.30) were true everywhere on the cut, we would im- 



mediately arrive at a paradox [ }40| . Consider the function f(ch(p) analytic 
in the chip plane with the cut from —1 to — oo. On the cut |/| 2 < and 
/(l) = 1. This is consistent with the maximum modulus theorem only at 

2ni < N c (2.31) 

what is not the case in our world. The more detailed analysis |Q shows that 
it is possible to obtain similar inequalities (with the constant smaller than 
2) using weight functions that are rather insensitive to the threshold region, 
and the paradox remains. 

HQET can also be used for description of heavy to light transitions ]41| 
and rare B decays ^] . In these cases the heavy quark spin symmetry is not 
so restrictive, and more form factors are necessary. Relations between B and 
D decays can be established using the heavy quark spin-flavour symmetry 
and the isospin symmetry. Inclusive heavy to light decays are considered 
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in pq |; they are described by several structure functions obeying sum rules 
in the deep inelastic region. 

An interesting approach in which the parameter ( y "^ ax ) 2 = ^ ™^™ c ) is 

considered small was proposed in f43| . To the leading order in this parameter, 
the quark-hadron duality is perfect: the b — > c W quark decay rate is equal 
to the B — > DW meson decay rate (see ( 2.26 )). This is true for all heavy 
quark polarizations; in particular, the hadronic tensor <B\J\XxX\ J + \B> 
summed over ground state mesons X = D, D* is equal to the corresponding 
quark tensor summed over c polarizations. 



3. Baryons 

For ground-state baryons, the light quark spins can add giving j p = + or 
1 + . In the first case their spin wave function is antisymmetric, the Fermi 
statistics and the antisymmetry in color require an antisymmetric flavour 
wave function. Hence the light quarks must be different; if they are u, d then 
their isospin I — 0. With the heavy quark spin switched off, we have the + 
1 = baryon Aq. If one of the light quarks is s, we obtain the isodoublet Sq 
forming together with Aq the SU(3) antitriplet. In the 1 + case the flavour 
wave function is symmetric; if the light quarks are u, d then their isospin 
1=1. So we have the 1 + isotriplet Eq; with one s quark — the isodoublet 
Sq; with two s quarks — the isosinglet £Iq. Together they form the SU(3) 
sextet. With the heavy quark spin switched on, the scalar baryons Aq, Sq 

become | + ; the vector baryons form degenerate | + doublets Eq, Eq; 
S Q' s q; n *Q- 

Baryons and their form factors in HQET were considered in p4|p5pq |. 
Two-point and three-point HQET sum rules were investigated in EW. 

Baryon currents with the scalar heavy quark have the form J s = 
£ ahc (q Ta CTTq b )Q c s where q T means q transposed and C is the charge con- 
jugation matrix (because q T C is transformed like q under the action of the 
Lorentz group). Here r is a flavour matrix, symmetric for + baryons and 
antisymmetric for 1 + ones. We shall abbreviate it to J s = (q T CTq)Q s . A 
light quark pair with j p = + corresponds to the current a = q T C n fc l q, and 
with 1 + — to a = q T C^jq (one can easily check it using the P-conjugation 
q — * 7o<z)- It is also possible to insert 70 into these currents without chang- 
ing their quantum numbers. So, the scalar heavy quark currents with the 
quantum numbers of Aq, Eq are Ja s = aQ s , J Ss = aQ s . 

With the real- world spin h heavy quark, the current J= aQ has the spin 
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the current j = aQ contains spin ~ and spin | components. The part 
J3/2 = J+ j77 ' J satisfies the condition 7- J 3 / 2 = and hence has the spin |. 
The other part J 1/2 = -§77 1 J= 3775 J1/2, J1/2 = « ■ 77sQ has the spin \. 
Finally we obtain the currents J= (q T CTq)T'Q with the quantum numbers 
ofAg, Eg, E Q Q 

J A = (<Z T C 75 g)g, JS = (g T C79)-7 75Q, %. = (q T Cjq)Q+^(q T Ciq)-jQ, 

(3.1) 

and similar currents with the extra 70 inside the brackets. 

x 





Figure 4: Correlator of two HQET baryonic currents 

The correlators of two baryon currents with the scalar heavy quark have 
the structure (Fig. |j) 

i<TJ As (x)JiM> = S(x)Ttt + tU a (x ), (3.2) 
i<TJ^(x)j+ SJ {0)> = SijSfflTlT+TlLxixo). 

From now on we shall for simplicity assume the normalization Tr r + r = 1. If 

we denote <0|Ja|Aq,0 + > = /a,o+j <0|j s |Sq, l + > = / s ,i+e, then the baryon 

contribution to pa.s(^) is f\ — ea,e)- 

Now let's switch the heavy quark spin on. The correlators are (Fig. ^) 

n=(r 1 i±^F' 2 W (3.3) 



where tensor indices may be contracted between IT S and T' 12 . The same 
relation holds for H(t), H(u>), and p(ui). For Aq, Eq, Eq we obtain 

l+7o 
Pa = — ^ — PAs, 
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1 



PT, 
PS* 



7*75 — 

Sw — 
l + 7o 



7o , „ 1 + 7o 

i IjlbOrjPSs = 3 — - — pss 



(3.4) 



g7i7i' 



l+7o 



>y + g7i7j ) Ps s - 



Si'j'PZs 



If we denote <0|ja|Aq, | > = / A i + u, <0|j s |Sq, ^ > = / E> i+U, 
<0|J s »|Eq,| + > = f s » 3 + u, then the baryon contributions to (3^) are 
^/^ + ^/ 2 , + ^-?e), ^ (6 l3 + | 7i7i ) 4 r *( w -e E ). 
Aq is a spin symmetry singlet, therefore there are no interesting predictions of 
the spin symmetry in this channel. Baryons in Sq and £q channels are degen- 



erate: e s i + 



£s,i+, and 75/2,1 + 



E.1+- 



Note that 



both sides of the definitions <0\j\B> = fsu get the same factor y2m when 
going from the relativistic normalization to the nonrelativistic one. Therefore 
the QCD quantities fs = Jb don't depend on m (compare with (2.9)). 




Figure 5: Correlator of two HQET baryonic currents and a heavy-heavy 
current 

The correlators of two baryon currents with the scalar heavy quark and 
the scalar heavy-heavy current J s have the structure (Fig. [^) 

i 2 <Tj As (x 2 )J s (0)j{ s (x 2 )> = Tr T+r 

OO 

dt 2 5{x 2 - v 2 t 2 ) J dt 1 5(x 1 +v 1 tl2)K A (t 2 ,t 1 ), (3.5) 

^ 2 <TJ^(x 2 )J s (0)^(x 2 )> = Trr+r 
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dt 2 5(x 2 - v 2 t 2 ) J dtiS(xi + Vitl2)Ksnu(h, h), 
o o 
Kznv = ^E||e2|| M ei||i, + Ke±5±h V , 

where ei|| = {v 2 — chipvi) / ship, egii = — (wi — ch(pv 2 )/ ship are the 
Eq polarization vectors in the scattering plane, Sj_^ v = "^2 e± fl e± l/ = 
[chip(v ltl v 2u + v 2fi v ll/ ) - v x ^v lv - V2ftV2v] / sh 2 ip - g^. 



According to the rules 26 , the transition Aq — > Aq is described by one 
form factor = j 2z = 0); Aq — > £q is forbidden by naturalness; £q — > £q 
is described by two form factors (ji z = j 2z = and ±1): 

<A Q \J S \A Q > = a(ch^), (3.6) 
<Sq|J s |Sq> = Cs^w/e^eiv, £ Smi/ = (ch y)e 2 | | M ei| |^ + £r;_L(cb.<p)£i_ /1 ,,. 

The contribution of Aq, Sq to Pa,e||,e_l(w2, ^i) is /a,e£a,e||,e_l£(w2 - 
£a,e)^(wi — £a,e)- The spectral density of the correlator of the currents jj^e^ 
with some specific polarizations ea.,2 is PT,^& 2 ^e-\ v ; the £q contribution to it 

is /ISs/u^ei!,. 

Now let's switch the heavy quark spin on. The correlators are (Fig. ||) 

K=(t>^T 1 -±^ 1 \k s (3.7) 



Let's introduce the currents ju — {qTC 1 q T )QB, B — V u. Rewriting (3.1) in 
the covariant form J s = (Q 1 C^^q)T'^Q, J^* v = (Q T Cj f _ l q)T'^Q, we obtain 
B\ = u, £>s M = — (jfj,+Vfj,)u, By,*^ = Ufi. The spectral density of the correlator 

is (j3 2 1+ 2 V2 r 1+ 2 V2 Bij p s (where tensor indices may be contracted between 

p s and Bi :2 ); the baryons' contribution to it is fB 2 fBi<B 2 \ J\Bi>S(u> 2 — 
£2)^(^1 - £i). Hence we obtain Qppl 



<Aq|J|Aq> = £ A (chip)u 2 r Ul , 

<Eq|J|Eq> = -^E /1I /W 2 (7 At + v 2fi )j 5 r(j„ - vi v )j 5 u 1: 
<Sq|J|Sq> = ■i^E/*wU2/*r(7 l , - vi v )^ 5 ui, (3.8) 

<Sq|J|Sq> = tepvUi/jFui,,. 

The result for Aq is particularly simple because light fields have j p = + , 
and the spin of Aq is carried by the heavy quark. 
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At the point ip = 



a(l) = &ll(l)=6a(l) = l 



(3.9) 



(at this point = Sij because there are no selected directions). 

Inclusive Aq decays were treated in With the scalar heavy quarks, 

the matrix element of the decay Aq — > X has the structure M — g(p2(X)ipi 
where phi\ is the scalar Aq wave function. The Aq decay rate is given 
by ( 2.22j ) with the structure function w(e, chip) = ^2 X (p2(X)cp2(X)5(ex — e) 
obeying the Bjorken sum rule (2.23). Transitions to the excited baryons and 
their contribution to the Bjorken sum rule were considered in 48,4^. 

Polarization effects in Aq decays were discussed in p3]. Heavy to light 

they are described by several 
s and sum rules for them in the 
a]. At a large number of colors, 
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baryon transitions were considered in 
form factors. Inclusive heavy to light deca 
deep inelastic region were investigated in J 



baryons are bound states of a chiral soliton and a meson; this model was 
considered in Kill ■ 



4. Renormalization 

Rcnormalization properties (anomalous dimensions etc.) of HQET are dif- 
ferent from that of QCD. The ultraviolet behavior of a HQET diagram is 
determined by the region of loop momenta much larger than all character- 
istic scales of the process but much less than the heavy quark mass which 
tends to infinity from the very beginning. It has nothing to do with the ul- 
traviolet behavior of the corresponding QCD diagram with the heavy quark 
line which is determined by the region of loop momenta much larger than the 
heavy quark mass. In the conventional QCD the first region produces hybrid 



logarithms |52 53 , and the problem of summation of these logarithmic cor- 
rections is highly nontrivial. In HQET hybrid logarithms become ultraviolet 
logarithmic divergences governed by the renormalization group with corre- 
sponding anomalous dimensions. For example, correlators of QCD meson 
currents contain large hybrid logarithmic corrections a s log — . Correlators of 
the corresponding HQET currents contain instead corrections a s log — where 
li is the normalization point. The dependence on fi is determined by the 
currents' anomalous dimensions; the corrections are small at /i ~ u). 



HQET is closely related to the theory of Wilson lines in QCD |54|. As 



follows from the Lagrangian ( |l.2[ ) , the static quark propagator in a gluon field 
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is the straight Wilson line 



S(x) = -itf(x )S(x)Pex.pig / A^dx^. (4.1) 



The effective Lagrangian identical to the HQET Lagrangian (1.5) (strictly 
speaking, with the scalar static quark Lagrangian ( 1.1 0| ) ) was proposed in 



for investigation of Wilson lines. Their renormalization properties were con- 
sidered in Q. 

One-loop renormalization of straight Wilson lines (static quark propaga- 
tors) and cusps on them (heavy-heavy velocity changing currents) are known 
from J54|. Two-loop calculatio n fl57| for straight Wilson lines is incorrect; the 
correct result was obtained in ]58|]. It was also obtained in [|9| starting from 
the on-shell renormalization of the QCD heavy quark propagator at finite m, 
and in |]60| , |6l| in the HQET framework. Two-loop renormalization of a cusp 
on a Wilson line was first considered in but the authors were unable to 
get rid of all double integrals. A simple result containing only simple single 
integrals was obtained in |^2| . The attempt in the HQET framework was 
unsuccessful: the result contains a double integral (with a variable undefined 
in the paper); some other integrals are in fact equal to or each other. 

One-loop renormalization of the heavy-light bilinear current in HQET was 
first considered in |5^j5^| . Two- loop corrections were obtained in [|60[|6l| (in 
the second paper, a different external momentum configuration was chosen 
which made calculations more difficult). Four-quark operators with two static 
quark fields were also investigated in [pl| . One-loop renormalization of baryon 
currents was considered in p7| . 

We shall use MS scheme, the space dimension D = 4 — 2e. The HQET 



lagrangian (1.2) expressed via bare fields and couplings is 



-fil^Gl^ + —{d^f + (ghost). (4.2) 



L = Q+i{d-ig b A a b t a ) Q b + q b i{d-ig b A a b t a )q 

The bare quantities are related to the renormalized ones as 

Q b =JI^Z Q /2 Q 1 q b =T+- £ Z 1 q / 2 q, ^=TT e A%, 

g b = JI e Zl /2 g, a b = Z A a, (4.3) 

where Z q , Z Al Z a are the same as in QCD with ni light flavours (there are 
no static quark loops), and /I 2 = /i 2 e 7 /47r, /i is the normalization point. 
The static quark field renormalization constant Zq is determined from the 
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requirement that the renormalized propagator S(ui) — Sb(uj)/ZQ is finite. 
If we denote the sum of bare one-particle^irreducible static quark diagrams 
— i£f,(w), then the propagator Sb{u) = So(u>) + So(u/)E;,(a;)So(w) + "" = 
l/(w-E 6 (w)). 

k 




^> =^ 



Figure 6: HQET mass operator 

The bare one-loop HQET mass operator (Fig. |^) in the Feynman gauge 
ab = is 

d D k 1 



^^^'W^TM (4 ' 4) 



where Cf — ~2N~' ^ variant of the Feynman parametrization 

oo 

1 T{a + 0) f yP^dy 



a a bP T(a)T((3) J (a + yb) a +P 



(4.5) 



is used to combine a square denominator a with a linear denominator b; the 
parameter y has the dimension of mass. We have 

oo 

d D k 1 

{2n) D (k 2 + yk + yujf 



E(w) - -iC F g 2 b I dy / ,1.2 , , „„/,2 ■ ( 4 ' 6 ) 



The denominator is equal to fc' 2 — \ + uiy, k' = k + \v. Using the standard 
formula 

d D k 1 _ t (-l)"r(n-f)(a 2 ) D / 2 -» 

(2tt) b (fc 2 - a 2 )« ~ (n- l)\(4n) D / 2 ' ( ^ 

we obtain 

OO j-y j 2 2 

^H-(§C r ( 2 -f)/(T-^) ^ 
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The integral 



y a (ay + b)f 3 dy 



6 a+/3+l r (_ x _ a _ + a J 



,«+l 



r(-/3) 



(4.9) 



is calculated using the substitution y = ~ (A — l). Finally. 



E(w) 



(4tt) d /2 



2(-2w)^- J r(3-z))r(f -l) 



= 1 + C f 



27re' 



(4.10) 



Requiring the finiteness of S(uS) — Sb{u)/ZQ with a minimal Zq = 1 + c— , 
we find 



(4.11) 




Figure 7: Heavy-light vertex 

Now we shall consider renormalization of a heavy-light bilinear current. 
The bare current % = Qb^Qb is related to the renormalized one as % = 
~p~ 2e ZjJ, or QTq = Zy J where Zj = Z^ 2 Zq 2 Zr- Then the matrix element 

r = <Q\QTq\q> = Zr<Q\j\q> where the matrix element of J is finite. The 
vertex T does not include corrections on the external legs because it con- 
tains renormalized fields. We shall calculate it in the one-loop approximation 
(Fig. 0) in the Feynman gauge. We are interested only in the ultraviolet di- 
vergence of the one-loop diagram that does not depend on external momenta. 
At zero external momenta we have 



r 



l-iC F glJ 



d u k fc 7o 



(2ir) D (k 2 ) 2 k 



= r 



1 - iC F gl 



d D k 1 
(2n) D (k 2 ) 2 



(4.12) 



because k = fco7o — fc • 7 and the integral with k vanishes due to the sym- 
metry. If we started from an infrared regularized matrix element (e. g. with 
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nonzero external momenta or gluon mass) we would obtain an integral with 
the same ultraviolet divergence but infrared safe. Separating the ultraviolet 
pole we have Zr = 1 + Cf-^t-. Using also (4.11) and the standard QCD 
renormalization constant Z q = 1 — Cf^; in the Fcynman gauge, we obtain 
the gauge invariant renormalization constant J52 53 4| 



1 



3C F ^. 

Site 



(4.13) 





Figure 8: Heavy- heavy vertex 

Similarly, the vertex A = <Q2\Q2Qi\Qi> m the one loop approximation 
(Fig. |^) in the Feynman gauge is 



f d D 

1 - iC F g 2 chip J — 



d u k 1 



D /j 2 (wifc + CJi)(w2fc + CJ 2 ) 

l 

2 . 



(4.14) 



1 — 2iCF9b ch</? J dx J dy 
o o 

1 



d D k 

(27T) 



D 



[k 2 + y(xv\ + (1 - x)v 2 )k + y{xu>i + (1 - x)uj 2 )] 3 



= 1 - 



1 OO 



■T(l + e)chtp J dx J dyy~ £ 



o o 

[\{x 2 + (1 - x) 2 + 2x(l - x) c\\ip)y - xu 1 - (1 - x)ui 2 ] 



-l-e 



(4tt) d / 2 V ' V ' Y J [x 2 + (l-x) 2 + 2x(l-x)ch<p] 1 - 

o 

Retaining only the ultraviolet 1/e pole and using the substitution x = i(l 
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zcth ^), we obtain 



+ th£ 



1 — Cf — — ethos 

2tT£ ^ 



1 + z 2 



(4.15) 



-th* 



and finally Za = 1 — Cf^^^cth^. Using ( 4.11 ) we obtain the gauge- 
invariant renormalization constant Zj = ZqZ^ [p3j53,p^,|l8f 



Zj = l- C F — ((pcthip - 1). 

27T£ 



(4.16) 



It is equal to 1 at <p = because then there is no cusp on the heavy quark 
world line. 

Knowledge of the renormalization constants Z allows us to determine the 
dependence of the renormalized operators on the normalization point /i using 
the renormalization group equations. Up to two loops, MS renormalization 
constants have the form 



47T £ V 47T / V £ 2 e J 

1 /2 

From gi — ~p £ Z a g = const we obtain the evolution of a s (/x) 



(4.17) 



I log a s 



-2( £ + /?(a s )),/?(a s ) = ±^ = 



that /?i 



^ c - In/, /3 2 



3 iv c 



1 

?2 



d log 

• • •. It is well known 
m. Substituting the 



form (4.17) for Z a we see that c\ — —fix, C22 = fit, C21 = — \fl2, i- e. C22 
is not independent. Similarly, the /Lt-dependence of any operator j is usually 

d log 2j 
d log /A 



characterized by its anomalous dimension 7^ 



Substituting the form (|4.17| ) for we see that ci = —571, C22 = 5(7? + 
2/3i7i), C21 = — j72, i- e. again C22 is determined by one-loop quantities. 
In the case of a non-gauge-invariant operator, Z also depends on the gauge 
parameter a(/x), and the formulae become more complicated. Solving the 
renormalization group equation we obtain 



-72 



4-tt ; 



1 + 


( 12 


71/^ 








4tt 



(4.18) 



where j is a renormalization group invariant. This formula allows us to re- 
late j(fJ-i) to j(p2)- Here we present for reference the two-loop anomalous 
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dimensions of heavy-light and heavy-heavy currents [16006 



7j 
7J 



4-7T 
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CfOc s 



((/9Cth(/3 — 1) 



- ni — {ipct\\ip - 1) 



Cfcx 2 s 



(4.19) 



/l /67 7T 2 \. 



cth (p tp ctYiipdip + cth 2 93 / — V') cth V'^V' 



ship f ip cth ip — 1 



-SI*) 



sh 99 - sh -0 sl1 ^ 



Until now we discussed the renormalization inside HQET. But usually we 
are interested in matrix elements of QCD operators (e. g. weak currents). 
Therefore we have to discuss the relation of QCD operators to their HQET 
analogues. Operators in HQET differ from those in QCD starting from the 
one-loop level even if written in the same form via the fields because their ma- 
trix elements are calculated using different Feynman rules. A QCD operator j 
matches the corresponding HQET operator AJ if they give identical physical 
(on-shell) matrix elements between states suitable for HQET treatment (with 
residual momenta much less than to). In order to calculate on-shell matrix 
elements we have to use the on-shell renormalization scheme in which prop- 
agators in the on-shell limit are free. For the "massless" fields q, Q the bare 
on-shell propagators get no corrections because loop integrals are no-scale 
(ultraviolet and infrared divergences cancel). Therefore the on-shell renor- 
malized fields coincide with the bare ones: q = 



-V2, 



-1/2, 



Note that although the expressions for the renormalization constants Z qi 



Z ( 



are the same as above, all divergences in them are infrared ones because these 
Z factors relate renormalized (ultraviolet-finite) fields. For the massive quark 



— 1/2 

field we have Q = Z Q 1 Q os , Zq = 1+Cf 



3L + 4) , where L = log : 



->Q *■ 1 4tt \ e ■ -/1 — o y 

The infrared divergence of the on-shell massive quark propagator Zq is the 
same as that of the static quark propagator Zq . 



For the heavy- light bilinear currents we have j = Z T Z, 



Qos r, 7o 



— Z r Zf_ 



-1/25-1/2 



<Q\j\q> = Z^Z ( 



<2os r <?o 
-1/2 r,-l/2 



Hence the on-shell matrix elements are 



r, <Q\j\q> 



-1 y-1/2 5-1/2 



z r z. 



r, where the 
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mv+k 




Figure 9: Light- light — > heavy-light matching 



proper vertices T, T are depicted on Fig. |[ We obtain the matching con- 
stant 



(4.20) 



Here ultraviolet divergences cancel in T/Zr, T / Zr by definition; infrared di- 
vergences cancel between these two expressions because the infrared behavior 
of QCD and HQET is identical; Aq is finite for the same reason. We choose 
all quark momenta in HQET to be zero; this corresponds to the heavy quark 
momentum mv (v — (1,0)) in QCD. Then the HQET loops (Fig. ^) vanish: 
r = 1. Let's calculate the QCD vertex T (Fig. ^) in the Feynman gauge 



r - iC F g 2 b 



r - 2iC F g£ / dx(l - x) 



dPk 7 M (fc + mv + m)rfc7 A1 
(2tt) d (fc 2 ) 2 (fc 2 + 2mvk) 



d D k 



(4.21) 



Jfj, (k' + m(l — x)v + m)r(fc' — mxv)j^ 



{k 12 — m 2 x 2 ) 3 
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where k' — k + mxv. The term with two k' in the numerator gives ^jyH 2 (D)T 
where 7^7^ = H(D)T; terms with one k' vanish. Terms without k! give 
±m 2 x7 /J (l — x + v)T^j^ = —m 2 x(2 ± H(D)x)T where the upper (lower) sign 
is for r anticommuting (commuting) with v and we have used the fact that v 
on the left may be replaced by 1 in the on-shell matrix element. Therefore all 
the terms have the common 7-matrix structure T; calculating the integrals, 
we have the vertex 



= 1 



a s / m" 

47T \ /i 2 



OL s I 171 

Of— —ft 



2 J dxx~ 2e 
'0 

1 



(l-x) 



4tt \n 2 J (l-e)(l-2e) 
The first divergence is ultraviolet: 



H 2 {D) 



4e 

H 2 {D) 1 - e 

4e e 



1 ± H{D) 
x 2 



± 



H(D) 



.(4.22) 



Z r = 1 + Cj 



47r 4e 



(4.23) 



By the way, the one-loop renormalization constant of the QCD bilinear quark 



currents is Z, 



ZqZ? 



Cj 



47r 4e 



the vector and axial current 



(H = ±2) anomalous dimension vanishes. Finally we obtain from ( 4.20| ) 
the matching [0 



QTq 



H 2 - 10 



L + -H 2 - HH 1 ± -H — 4 



where H' 



dH 



4tt V 4 " ' 4 2 

(4.24) 

, ;) . This equation holds separately for QCD currents with 
r (anti-) commuting with 70. If it does not have this property, it can be 
split into a commuting and an anticommuting part; it then maps onto a 
combinati on of two HQET currents. The logarithmic part of the matching 
constant ( 4.24 ) is determined by the difference of anomalous dimensions of 
the QCD and HQET currents; the non-logarithmic part should be included 
accoun t onl y when the two-loop anomalous dimension is also taken into ac- 
count |yj). _ 

Now we shall consider the current Q 2 TQi in the effective theory where Qi 
is a heavy quark with the mass to. We can go to the second effective theory 
in which this quark is also considered static. The matching is obtained by 
comparing the diagrams in Fig. 10. All external momenta are zero in the 



second theory, therefore the loops vanish. In the first theory, the heavy quark 
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Figure 10: Heavy-light — ► heavy-heavy matching 



has the momentum mv\ . The vertex is the matrix T times 
1 - iC F g 2 b 



d D k (k + mv\ + m)v2 
'"' 1 (2w) D k 2 (k 2 + 2mv 1 k)v 2 k 



(4.25) 



1 - 2iC F gi dy dx 



o o 



d D k —y/2 + m(l — x)v 1 V2 + mv2 
(2ir) D (k 12 — y 2 /A — m 2 x 2 — mxychf) 3. 



Now we calculate the loop integrals. The parametric integrals factorize after 
using the substitution y — 2mxz, and the x integrals are trivial. The upper 
(lower) sign is for T anticommuting (commuting) with v 2 . 



a 8 I m 



2\ -e 



? 1 _ 2£ xz-2(l-x)ch( P ± a; 

i+2cf 47t Uv y y (i+z 2 +2zch^^ 





(4.26) 



1 + 2C F - 



(z + 2chy± !)/(! -2 £ )-ch^/e 
(l + z 2 + 2zch( / 9) 1 + e 



The first ultraviolet divergent integral can be calculated by splitting the in- 
tegration region at a large A and ignoring e in the first region and 1/z in the 
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second one: 



z dz 



(1 + z 2 + 2zchip) 1+£ 



z dz 



1 + z 2 + 2z ch <p 



1 



dz 

yl+2e 



(4.27) 



log A — <p cth ip H log A. 

2e 



The second integral is convergent; we need it up to 0(e) because it is multi- 
plied by the infrared pole 1/e: 



oo 



dz 



o 

where 



(1 + z 2 + 2zchp) 1+e shp 



F(x) 



log(l + y) 



dy 



(4.28) 
(4.29) 



is the Spence function. Two Spencc functions in (4.28) are not independent: 
F(e 2v - 1) + F(e- 2 ^ - 1) = 2ip 2 . The vertex is 



1 + C *T. 



1 2 ip 

ip cth p — L + 2Lip cth ip + 2 + 2 ip cth ip ± 2— — 

e e ' ' shp 



- cth <p (F(e 2v - 1) - F(e- 2¥; - 1)) 



(4.30) 



The first divergence is ultraviolet: Zr = l + C^fj-; this is the renormaliza- 
tion constant of the hea vy-light current, and it in deed agrees with what we 
have found before ( 4.13 ). In the denominator of ( 4.20| ), we should use the 



renormalization constant of the heavy-heavy current -Za found before (4.16) 



The infrared divergence cancels as it should do, and we obtain the match- 

ing (Tig 



Q 2 TQ 1 = 



1 + C F — [ 2Lp cth p> +-L + 2p cth p ± 2-^- 
47r \ 2 shp 



+ cthtp (F(e 2v - 1) - F{e- 2ip - 1)) + 



Q 2 TQ 1 (4.31) 
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m x v x 



X™2^2 



k 




Figure 1 1 : Light-light — > heavy-heavy matching 



Finally we shall consider the current (3 2 r<3i where Qi,2 are the heavy 
quarks with the masses mi^. We can go to the effective theory in which both 
of them are considered static. The matching is obtained by comparing the 
diagrams in Fig. 11. All external momenta are zero in the effective theory, 
therefore the loops vanish. In QCD, the heavy quarks have the momenta 
mi 2^1 2 ■ The vertex is 



r 



d D k 7 M (fc + m 2 V2 + m 2 )r(fc + m\V\ + mi)^^ 



Y — HCf% / dx\dx 



k 2 (k 2 + 2m\V\k)(k 2 + 2ni2V2k) 
d D k 



(4.32) 



(2n)V 

7 M (fc' + m 2 {\ - x 2 )v2 - mxxxvx + m 2 )r(fc' + m^l - x^vi - m 2 V2 + mi)^ 

(fc' 2 - a 2 ) 3 ' 

where a 2 = m\x\ + m 2 x 2 + 2mim2X\X2ch. tp. Now we calculate the loop 
integrals. The parametric integrals factorize after the substitution x\ t 2 = 
a;(l±z)/2, and the x integrals are trivial (a 2 = mim2X 2 a+a-, a± = ch ^ y + 
zsh^, where = log^L): 



dz 



H 2 (D) 
8e(l -e, 



ch(/3 



e(l — 2e)a + a_ 
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+ 



chip + zshiP H{D)(l-z 2 ) 



(1 — 2e)a + a_ 16(1 — e)a + a_ 



(4.33) 



dz 



[l + z+\He-^(l-zf 



^ ^ a s H 



l-z< 



dz. 



At mi = (ip — > — oo) this vertex reduces to (4.22), and at m2 = (V 1 — > 
+oo) — to the mirror symmetric expression. In order to check this, we s hould 
take the limit before e — > 0. Ultraviolet divergences are the same ( 4.23 ), but 
the structure of infrared 1/e poles is different at zero and nonzero masses. 
The integrals in z are easily calculable: 



+i 



dz 



-l 
+ i 



dz 



(a+cj- 




(4.34) 
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+i 







sh^±£ 




-(¥>- 




2ih log -1— 


)■ 
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sh ^ — ip sh </5 






shy 


ch ip — ch V 1 








4 


4^ sh^ 




ch 


ip — chip (chip — chip) 2 



2(y5 sh + sh ip 
sh <p (ch ip — chip) 2 



The Spence functions here are not independent: F ^ ll^tj? J + ~^ — ¥") = 

| l°g 2 e^-e-f • Matching is determined by the formula similar to (4.20) but 
with AqiAq2', infrared divergences cancel as they should do, and we finally 
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obtain |E3U)4 



Q 2 rQi = AQ 2 TQ 1 + A 1 Q 2 v 1 TQ 1 + A 2 Q 2 Tv 2 Qi + A X2 Q 2 v 1 Tv 2 Q x , (4.35) 



A = 1 + C F ^ 

4-7T 



H 2 



— — h 2<p cth (/? — 3 ) L — 



H 2 



-HH' 



H ( ip ch^chi/i — 1 



4 

ipship 



ipship — ipship 
chip — chip 

1 



ship (chip — chip) 2 (chip — chip) 2 chip — chip 



-cth>p F 



(ip 



-,) + 2^1og^ 



^ ip chip chip — 2 ch 99 + 1 ^ V'shV' 



" Cf 4^ 



ch 93 — ch ip 
1 



ch y> — ch ip 



ip ( ch ip ch ip — 1 ^ ip ship — ip sh ip 



4 \ ship \ (chip — ch-0) 2 



ip ship 



(chip — chip) 2 chip — chip 



sh ip 



ch — ch -0 

ip sh ip — ip sh <p 
chip — ch ip 



1 



A 2 = -C F -± 

47T 

ip ship 



ip ( chip chip — 1 ipship! — ipship 
chip — chip 



ship \ (chip — chip) 2 



1 



(chip — chip) 2 chip — chip 



sh ip 



1 



A12 — Cf-t-—t 



ip ch ip ch ip — 1 



ipship — ipship 
chip — chip 

ipship 1 



Aw 4 \sh(/? (chy) — ch-0) 2 (chip — chip) 2 chip— chip J' 



where L = log 



The one- loop matching of the baryonic currents was considered in |4jJ . 
Now we are in a position to make some statements of Sec. ||, || more 
precise. The QCD meson constants (2.9) are related to the HQET constant 
by the matching (4.24): 



2/(m) 



a s (m) 



(4.36) 



where c = | for vector mesons and c — | for pseudoscalar mesons (if a fully 
anticommuting 75 is used). The HQET constant depends on the normaliza- 
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tion point as ( [4.18| - [PS| ): 



/( M ) = 7a7 2 ^(M)(l 



k ~ 12 



285 - 7?H 



107 



27/3i 

(4.37) 

There are two approaches to the b — > c weak decays: one-step match- 
ing and two-step matching |p!s| , |S5|j6^| , We are interested in hadronic 
matrix elements of the vector and axial weak currents j = cTb, T = 7^ 
or 7^75. These currents are defined in QCD at a high normalization point 
fi ~ mw- In the first approach, we use QCD at the scales from mw down 
to some not exactly definable border m ~ nib ~ m c - By a chance, the 
QCD anomalous dimensions of these currents vanish, and j(fn) = j(mw)- 
At n = m we perform the matching (4.35) to the HQET in which both b and 
c quarks are considered static. The vector current becomes a combination of 
cTb with r = 7^, Wfcp, and v cfJ .', the axial current — of the similar currents with 
the extra 75. Then we scale down to a typical hadronic /i using the HQET 
heavy- heavy anomalous dimension (4.1E). At this point we use the heavy 
quark spin symmetry, and express the matrix elements via the Isgur-Wise 
form factors. 

In the two-step approach, we use QCD from /i = mw down to /1 = m^. 



At this point we perform the matching (4.24) to the HQET-1 in which b is 
static while c is still dynamic. The vector current becomes a combination of 
cTb with r = 7 M and Vbfi (and the extra 75 in the axial case) . Then we use the 
HQET heavy-light anomalous dimension (4. IE) to scale these currents down 
to fx = m c . At this point we perform the matching (4.31) to the HQET-2 in 
which both b and c are static. We obtain a combination of cTb with r = 7^, 
Vb/j,, and v cfi (with the extra 75 in the axial case). Then we use the HQET 
heavy- heavy anomalous dimension ( 4.19| ) and the spin symmetry as before. 

In the one-step approach, we can't sum the a s log ^ corrections; we can 
do it in the two-step approach (even in the subleading order). On the other 
hand, the first matching in the two-step method gives a series in 2** because 
m c is the largest mass scale in the intermediate HQET. In the above descrip- 
tion all ^£ corrections were discarded, and this is not a good approximation 
in the real world. The first ^ correction can be included |^| (the leading 
a s log ^ corrections are summed in this term using the one-loop anomalous 
dimensions), but incorporating the second term would require a large work. 

The one-step matching seems more adequate in our world in which ^ is 
not very small and log ^ is not too large. It is possible to obtain the optimal 
combination of both approaches [0 . We expand the result of the one-step 
matching in 2U, Then we extract the zeroth term from the series, and replace 
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it by the result of the two-step matching. Wc also extract the first term, and 
replace it by the first power correction from the two-step matching. The 



errors of this procedure are of the order of a 2 s , or ^a s , or ( ^ ) a s log 



they all are small. 

References 

[1] H. Georgi. Proc. Theoretical Advanced Study Institute, ed. R. K. Ellis, 
C. T. Hill, J. D. Lykken, World Scientific (1992), p. 589. 

[2] B. Grinstein. High Energy Phenomenology, Proc. Workshop, Mexico 
City, cd. R. Huerta, M. A. Perez, World Scientific (1992); Annual 
Review of Nuclear and Particle Science 42 (1992). 

[3] M. B. Wise. Preprint CALT-68-1721, Pasadena (1991); N. Isgur, 
M. B. Wise. Preprint CEBAF-TH-92-10, Newport News (1992). 

[4] E. Eichten, B. Hill. Phys. Lett. B234 (1990) 511; B243 (1990) 427. 

[5] B. Grinstein. Nucl. Phys. B339 (1990) 253. 

[6] H. Georgi. Phys. Lett. B240 (1990) 447. 

[7] A. F. Falk, M. Ncubcrt, M. Luke. Preprint SLAC-PUB-5771, Stanford 
(1992). 

[8] M. Luke, A. V. Manohar. Phys. Lett. B286 (1992) 348. 

[9] M. J. Dugan, M. Golden, B. Grinstein. Phys. Lett. B282 (1992) 142. 

[10] N. Isgur, M. B. Wise. Phys. Lett. B232 (1989) 113; B237 (1990) 527. 

[11] H. Georgi, M. B. Wise. Phys. Lett. B243 (1990) 279; C. D. Caronc. 
Phys. Lett. B253 (1991) 408. 

[12] M. J. Savage, M. B. Wise. Phys. Lett. B248 (1990) 177; M. J. White, 
M. J. Savage. Phys. Lett. B271 (1991) 410. 

[13] E. Eichten. Nucl. Phys. B (19() Proc. Supplement) 4 (1988) 170; 20 
(1991) 475. 

[14] Ph. Boucaud, C. L. Lin, O. Pcne. Phys. Rev. D40 (1989) 1529; 
Erratum: D41 (1990) 3541; Ph. Boucaud, J. P. Leroy, J. Micheli, 
O. Pcne, G. C. Rossi. Preprint CERN-TH-6599/92 (1992). 




38 



[15] O. F. Hernandez, B. R. Hill. Phys. Lett. B237 (1990) 95; E. Eichten, 
B. Hill. Phys. Lett. B240 (1990) 193; J. M. Flynn, B. R. Hill. Phys. 
Lett. B264 (1991) 173. 

[16] L. Maiani, G. Martinelli, C. T. Sachrajda. Nucl. Phys. B368 (1992) 
281. 

[17] C. R. Alton, C. T. Sachrajda, V. Lubitz, L. Maiani, G. Martinelli. 
Nucl. Phys. B349 (1991) 598; C. Alexandrou, F. Jegerlehner, 
S. Giisken, K. Schilling, R. Sommcr. Phys. Lett. B256 (1991) 60. 

[18] A. F. Falk, H. Gcorgi, B. Grinstcin, M. B. Wise. Nucl. Phys. B343 
(1990) 1. 

[19] F. Hussain, J. G. Korner, K. Schilcher, G. Thompson, Y. L. Wu. Phys. 
Lett. B249 (1990) 295. 

[20] M. Ncubcrt, V. Rieckert, B. Stech, Q. P. Xu. Heavy Flavours, ed. 
A. J. Buras, M. Lindner, World Scientific (1992). 

[21] T. Manncl, W. Roberts, Z. Ryzak. Phys. Lett. B254 (1991) 274; 

M. Neubert. Phys. Lett. B264 (1991) 455; P. Ball. Phys. Lett. B281 
(1992) 133; G. Burdman. Phys. Lett. B284 (1992) 133. 

[22] T. Mannel, Z. Ryzak. Phys. Lett. B247 (1990) 412. 

[23] A. F. Falk, B. Grinstcin. Phys. Lett. B249 (1990) 314. 

[24] N. Isgur, M. B. Wise. Phys. Rev. D43 (1991) 819. 

[25] S. Balk, J. G. Korner, G. Thompson, F. Hussein. Preprint IC/91/397, 
Trieste (1991). 

[26] H. D. Politzer. Phys. Lett. B250 (1990) 128. 

[27] A. F. Falk. Nucl. Phys. B378 (1992) 79. 

[28] T. Mannel, W. Roberts. Preprint CEBAF-TH-92-21, Newport News 
(1992). 

[29] E. V. Shuryak. Nucl. Phys. B198 (1982) 83. 

[30] D. J. Broadhurst, A. G. Grozin. Phys. Lett. B274 (1992) 421; 

E. Bagan, P. Ball, V. M. Braun, H. G. Dosh. Phys. Lett. B278 (1992) 
457. 



39 



[31] A. V. Radyushkin. Phys. Lett. B271 (1991) 218; M. Neubert. Phys. 
Rev. D45 (1992) 2451; B. Blok, M. A. Shifman. Preprint 
TPI-MINN-92-32/T, Minnesota (1992); E. Bagan, P. Ball, 
P. Gosdzinsky. Preprint HD-THEP-92-40, Heidelberg (1992); 
M. Neubert. Preprint SLAC-PUB-5992, Stanford (1992). 

[32] B. Grinstcin, P. F. Mcndc. Phys. Rev. Lett. 69 (1992) 1018; 

M. Burkardt, E. S. Swanson. Preprint MIT CTP 2096, Massachusetts 
(1992). 

[33] T. Altomari, L. Wolfcnstcin. Phys. Rev. D37 (1988) 681; N. Isgur, 
D. Scora, B. Grinstein, M. B. Wise. Phys. Rev. D39 (1989) 799. 

[34] J. D. Bjorkcn. Preprint SLAC-PUB-5278, Stanford (1990); Proc. 18 
SLAC Summer Institute, SLAC-REPORT-378, Stanford (1991), p. 167. 

[35] J. D. Bjorkcn, I. Dunictz, J. Taron. Nucl. Phys. B371 (1992) 111. 

[36] M. Burkardt. Phys. Rev. D46 (1992) R1924. 

[37] J. Chay, H. Gcorgi, B. Grinstcin. Phys. Lett. B247 (1990) 399. 

[38] M. B. Voloshin. Preprint TPI-MINN-92/25-T, Minnesota (1992). 

[39] E. de Rafael, J. Taron. Phys. Lett. B282 (1992) 215. 

[40] D. J. Broadhurst. Preprint OUT-4102-38, Milton Keynes (1992), to be 
published in the Proc. Sotchi Conf. on High Energy Physics and 
Quantum Field Theory (1992). 

[41] N. Isgur, M. B. Wise. Phys. Rev. D42 (1990) 2388; G. Kramer, 
T. Manncl, G. A. Schuler. Zcit. Phys. C51 (1991) 649. 

[42] A. Ali, T. Manncl. Phys. Lett. B264 (1991) 447; G. Burdman, 
J. F. Donoghuc. Phys. Lett. B270 (1991) 55. 

[43] M. B. Voloshin, M. A. Shifman. Sov. J. Nucl. Phys. 47 (1988) 511; 
S. Rudaz, M. B. Voloshin. Phys. Lett. B252 (1990) 443. 

[44] N. Isgur, M. B. Wise. Nucl. Phys. B348 (1991) 276. 

[45] H. Gcorgi. Nucl. Phys. B348 (1991) 293; T. Manncl, W. Roberts, 
Z. Ryzak. Nucl. Phys. B355 (1991) 38. 



40 



[46] F. Hussain, J. G. Korner, M. Kramer, G. Thompson. Zeit. Phys. C51 

(1991) 321; F. Hussain, D. Liu, M. Kramer, J. G. Korner, S. Tawfiq. 
Nucl. Phys. B370 (1992) 259. 

[47] A. G. Grozin, O. I. Yakovlcv. Phys. Lett. B285 (1992) 254, B291 

(1992) 441. 

[48] N. Isgur, M. B. Wise, M. Youssefmir. Phys. Lett. B254 (1991) 215. 

[49] W. Roberts. Preprint CEBAF-TH-92-04, Newport News (1992). 

[50] T. Manncl, W. Roberts, Z. Ryzak. Phys. Lett. B255 (1991) 593; 
J. G. Korner, M. Kramer. Phys. Lett. B275 (1992) 495. 

[51] E. Jenkins, A. V. Manohar, M. B. Wise. Preprints CALT-68-1783, 
CALT-68-1809, Pasadena (1992); Z. Guralnik, M. Luke, 
A. V. Manohar. Preprint UCSD/PTH-92-24, La Jolla (1992); 
E. Jenkins, A. V. Manohar. Preprint UCSD/PTH-92-26, La Jolla 
(1992). 

[52] M. B. Voloshin, M. A. Shifman. Sov. J. Nucl. Phys. 45 (1987) 292. 

[53] H. D. Politzer, M. B. Wise. Phys. Lett. B206 (1988) 681; B208 (1988) 
504. 

[54] A. M. Polyakov. Nucl. Phys. B164 (1980) 171. 

[55] J. L. Gervais, A. Neveu. Nucl. Phys. B163 (1980) 189; I. Ya. Aref'eva. 
Phys. Lett. B93 (1980) 347. 

[56] V. S. Dotscnko, S. N. Vergeles. Nucl. Phys. B169 (1980) 527; 

R. A. Brandt, F. Neri, M. Sato. Phys. Rev. D24 (1981) 879; H. Dorn. 
Fortschr. Phys. 34 (1986) 11. 

[57] S. Aoyama. Nucl. Phys. B194 (1982) 513. 

[58] D. Knauss, K. Scharnhorst. Annalcn der Physik 41 (1984) 331. 

[59] D. J. Broadhurst, N. Gray, K. Schilcher. Zeit. Phys. C52 (1991) 111. 

[60] X. Ji, M. J. Musolf. Phys. Lett. B257 (1991) 409; D. J. Broadhurst, 
A. G. Grozin. Phys. Lett. B267 (1991) 105. 

[61] V. Gimenez. Nucl. Phys. B375 (1992) 582; Preprint INFN 873, Roma 
(1992). 



41 



[62] G. P. Korchemsky, A. V. Radyushkin. Nucl. Phys. B283 (1987) 342; 
Phys. Lett. B279 (1992) 359; G. P. Korchemsky. Mod. Phys. Lett. A4 
(1989) 1257; W. Kilian, P. Manakos, T. Mannel. Preprint IKDA-92/9, 
Darmstadt (1992). 

[63] X. Ji. Phys. Lett. B264 (1991) 193. 

[64] M. Ncubcrt. Nucl. Phys. B371 (1992) 149; Phys. Rev. D46 (1992) 
2212. 

[65] A. F. Falk, B. Grinstcin. Phys. Lett. B247 (1990) 406. 



42 



